Level-Spacing Distributions of Planar Quasiperiodic Tight-Binding Models 
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We study the statistical properties of energy spectra of two-dimensional quasiperiodic tight-binding 
models. We demonstrate that the nearest-neighbor level spacing distributions of these non-random 
systems are well described by random matrix theory. Properly taking into account the symmetries 
of models defined on various finite approximants of quasiperiodic tilings, we find that the underly- 
ing universal level-spacing distribution is given by the Wigner-Dyson distribution of the Gaussian 
orthogonal random matrix ensemble. Our data allow us to see the differences between the Wigner 
surmise and the exact level-spacing distribution. In particular, our result differs from the criti- 
cal level-spacing distribution computed at the metal-insulator transition in the three-dimensional 
Anderson model of disorder. 
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Following the pioneering works of Wigner and Dyson 
random matrix theory (RMT) has been successfully 
applied to investigate a great variety of complex systems 
such as nuclear spectra, large atoms, mesoscopic solids, 
and chaotic quantum billiards In such systems, it 

has been shown that spectral fluctuations can be modeled 
by universal level-spacing distributions (LSD) such as, 
e.g., -Pgoe(s) for the Gaussian orthogonal random matrix 
ensemble (GOE) §. 

A very natural application of RMT concerns disordered 
systems Q . It has been shown that the metal- insulator 
transition (MIT) in the three-dimensional (3D) Ander- 
son model of localization is accompanied by a transi- 
tion of the LSD P(s) Here, s denotes the energy 
spacing in units of the mean level spacing A. In the 
metallic regime, P(s) closely follows the Wigner surmise 
Pw(s) = 7rsexp(— 7rs 2 /4)/2, which is a good approxi- 
mation of Pgoe(s). On the insulating side, P(s) is given 
by Poisson's law Pp(s) = exp(— s). One important differ- 
ence between the two distributions is their small-s behav- 
ior: Pw(s — > 0) ~ ns/2 and Pp(s — ► 0) ~ 1, indicating 
level repulsion and clustering, respectively. At the MIT, 
where the eigenstates are multifractal |l(|, another LSD, 
P c (s), has been observed 

Multifractal eigenstates — neither extended nor ex- 
ponentially localized — have also been found in tight- 
binding (TB) models of quasicrystals. In fact, these seem 
predominant in ID and 2D jll); in 3D, the attainable sys- 
tem sizes are yet too small for definite statements [ p.2| . 
The multifractality is assumed to be connected to the 
unusual transport properties of quasicrystals |fL3|| , e.g., 
their resistivity increases considerably with decreasing 
temperature and improving structural quality of the sam- 
ple. Thus, one may speculate that the LSD in quasiperi- 
odic models is also distinct from the Wigner and Poisson 
behavior. 

Quasicrystals lack the translational symmetry of peri- 
odic crystals, but still retain long-range (orientational) 
order and show non-crystallographic symmetries incom- 



patible with lattice periodicity. Thus, they constitute 
a class of materials somewhere in between perfect crys- 
tals and amorphous systems. Besides quasicrystals with 
icosahedral symmetry [flif , which are aperiodic in any di- 
rection of the 3D space, also dodecagonal fl5|| , decagonal 
ipjSf , and octagonal phases have been found, which 
can be viewed as periodic stacks of quasiperiodic planes 
with 12-, 10-, and 8-fold symmetry, respectively. Struc- 
ture models of quasicrystals are based on quasiperiodic 
tilings which can be constructed, e.g., by projection from 
higher-dimensional periodic lattices |lq| . We emphasize 
that such quasiperiodic tilings, albeit yielding perfect ro- 
tationally symmetric diffraction patterns, exhibit n-fold 
rotational symmetry in a generalized sense only. In par- 
ticular, there need not be a point with respect to which 
the tiling has an exact global n-fold rotational symmetry. 
If such a point exists, it is unique. 

In order to understand the transport properties of 
quasicrystals JDfl , TB models defined on such aperi- 
odic tilings (notably the Penrose tiling) have received 
considerable attention Jll],[l9| For a TB model de- 
fined on the octagonal (Ammann-Beenker) tiling |2q] , the 
LSD has also been used to classify the spectrum [ES E5[ . 
For periodic approximants, level repulsion was observed 
p3| , p4j , and P{s) was argued to follow a log-normal dis- 
tribution p4| . However, a calculation for finite patches 
with an exact 8-fold symmetry yielded level clustering 
II- 

In this Letter, we show that these somewhat diverg- 
ing results become comprehensible when one realizes that 
the tilings of Refs. |^,^| still retain non-trivial symme- 
tries, such as a reflection symmetry for the standard pe- 
riodic approximants. In order to obtain the underlying 
universal LSD, one should consider the irreducible sub- 
spectra separately, or break the symmetry by, e.g., cither 
choosing patches without symmetry, or imposing suitable 
boundary shapes as in quantum billiards, or introducing 
disorder. Properly taking this into account, we find that 
the underlying LSD of these non-random Hamiltonians 
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is neither Pp(s) |p5| . nor log-normal pil], nor P c (s) but 
rather -Pgoe(s)- The accuracy of our data further allows 
us to show that it is also not Pw(s) although this is as 
usual a reasonable approximation J2| . We emphasize that 
our results apply to all planar tilings mentioned above. 

Let us reconsider [^3j-||] the octagonal tiling consist- 
ing of squares and rhombi with equal edge lengths as 
in Fig. |l](a). Besides the projection method mentioned 
above, one may also use the self-similarity of the tiling 
to construct ever larger patches by successive inflation 
steps ^7j. E.g., the patch in Fig. |l|(a) corresponds to 
two inflation steps of the inner shaded octagon. On this 
tiling, we define the Hamiltonian H = j) K)0'l with 
free boundary conditions, \i) denotes the Wannier state 
at vertex z, and runs over all pairs of vertices con- 
nected by an edge of unit length. 

We diagonalize the Hamiltonian by standard methods, 
and study the LSD of the full spectrum. Due to the bi- 
partiteness of the tiling, the energy spectrum is symmet- 
ric about E = 0. Furthermore, it has been shown that 
a finite fraction of the states are degenerate at E = 
jl9],^2 - 24 1 . These correspond to confined states Jl9| lim- 
ited to certain local environments, do not contribute to 
the LSD, and we neglect them. In agreement with previ- 
ous calculations |^3| , we find that the integrated density 
of states (IDOS) is very smooth. This is different from 
ID quasiperiodic systems which typically have singular 
continuous spectra fu] . Nevertheless, the IDOS is not 
strictly linear as required by RMT, so we "unfold" the 
spectrum by fitting the IDOS to a cubic spline Q and 
use Si = N av (E i+ i) — N av (Ei) for the level-spacing at 
the ith level with A^ av the smoothed IDOS. We remark 
that the unfolded LSD is not a bulk quantity since A -1 
is proportional to the system size. In what follows, we 
shall always consider instead of P(s) the integrated LSD 
I(s) = f P(t)dt which is numerically more stable §,§. 

Fig. H(a) shows I(s) obtained for a large octagonal 
patch with 157369 vertices corresponding to three more 
inflation steps of Fig. |l|(a). At first glance, I(s) seems to 
be close to the integrated Poisson law Ip(s) = Pp(s), as 
observed in Ref. P5fl. However, this patch has the full 
Ds-symmetry of the regular octagon, hence the Hamil- 
tonian matrix splits into ten blocks according to the ir- 
reducible representations of the dihedral group D%: us- 
ing the rotational symmetry, one obtains eight blocks, 
two of which split further under reflection, while the 
remaining six form three pairs with identical spectra. 
This gives a total of seven independent subspectra. As 
with the confined states, we neglected the exact degen- 
eracies induced by symmetry in Fig. |^(a), since they 
only contribute to P(0). The integrated LSD of the 
seven independent subspectra are shown in Fig. |^(b). 
Each subspectrum matches the integrated Wigner sur- 
mise Iw(s) = exp(— 7rs 2 /4) to very good accuracy. From 
RMT it is known M that the LSD of a superposition 



of k independent spectra, each of which obeys Wigner 
statistics, is given by P^ (s) = ^ [erfc(^|)] fc with 
erfc(i) the complementary error function. For large k, 
P^\s) approaches the Poisson law Pp(s). In Fig. 0(a), 
we therefore also included the integrated LSD iff (s) of 
k = 7 Wigner spectra. The data clearly fits this curve 
better than Ip(s). This explains why in a previous cal- 
culation [^U a Poisson-like distribution was found. We 
also note that our data do not follow the integrated LSD 
I c (s) found at the Anderson MITJ§. 

Upon closer inspection of Fig. g(b) , we see that there 
are only very small differences between the seven inte- 
grated LSD, whereas there are slightly larger deviations 
to Iw( s )- In Fig. ||, we show the small- and large- s be- 
havior in more detail, restricting ourselves to data from 
one irreducible sector. We include data for patches of 
different sizes, corresponding to two, three, four, and five 
inflation steps of the inner shaded octagon of Fig. |l](a) 
with 833, 4713, 27137, and 157369 vertices, respectively. 
The convergence to /w(s) with increasing patch size is 
apparent both for small and large s in Fig. ^. However, 
the above-mentioned small deviations from iw( s ) still 
persist, even though the finite-size dependence is already 
very small. Therefore, we also consider in Fig. ^| the ex- 
act integrated LSD Igoe(s) 0. Although the Wigner 
surmise is usually a sufficient approximation of the exact 
LSD, we show in the inset of Fig. ^ that our data follow 
the exact curve /goe(s). Thus, the small deviations seen 
in Fig. ||, are due to the differences between Iw(s) and 
^goe(s). 

We can also approximate the octagonal tiling by 
patches which do not have exact symmetries. For in- 
stance, in Fig. |l](b) we show such a square-shaped patch 
cut out of the octagonal tiling. Although the quasiperi- 
odic eightfold order is restored in the infinite patch, there 
is never any exact symmetry present in the finite ap- 
proximants. The LSD is of the Wigner-type as shown in 
Fig. ||(c) for patches with side lengths L = 40, 60, and 
80, corresponding to 1980, 4392, and 7785 vertices, re- 
spectively. Thus, contrary to the case of a simple square 
lattice exhibiting level clustering, we find level repulsion. 
Again, we observe small deviations from I-w(s) which can 
be explained as previously when using /goe(s)- If one 
uses square-shaped approximants with symmetries, for 
instance squares centered around the eightfold point of 
the patch in Fig. ^(a) or the standard periodic approx- 
imants used in Rcfs. |2^j24|], the LSD is again a super- 
position of the LSD of the irreducible subspectra. Thus, 
approaching the infinite tiling by square-shaped patches 
only slightly shifted with respect to each other may give 
quite different LSD. We have obtained the same results 
also for circular patches. In this case, one can have ei- 
ther Dg, or reflection, or no symmetry, depending on the 
choice of the center. Thus, the LSD is well approximated 
by p w( s )> or p w ( s )' or p w( s ) = Av(s), respectively. 



2 



A different way of excluding any symmetries is given by 
choosing patches with special boundary shapes. In fact, 
this is well known in the context of quantum billiards 
where it has been used to construct quantum chaotic 
motion pLp|. One of the most prominent examples is 
the Sinai billiard |||5|, which consists of 1/8 of a square 
and a circular arc centered in the midpoint of the square. 
Due to these boundary conditions, the LSD follows the 
Wigner surmise even for free electrons [Q instead of a 
Poisson law which is found for integrable motion in sim- 
ple square or circular billiards In Fig. |l](c), we show 
a Sinai billiard-shaped cut of the octagonal tiling. Mov- 
ing Sinai's billiard table across the octagonal tiling, we 
can now generate many different patches. However, in 
contrast to the square-shaped boundary, we never find 
a case that retains any of the Ds-symmetries. We com- 
puted I(s) for quasiperiodic billiards with L — 70, 80, 90, 
100, and 110, corresponding to patches with 2416, 3146, 
3969, 4892, and 5905 vertices, respectively. The results 
presented in Fig. ||(d) follow I\y(s), and, again, are even 
closer to /goe(s). 

We emphasize that, apart from statistical fluctuations 
at small and large values of s as shown in Fig. |^, there 
is no systematic size-dependence of I(s). This is in con- 
trast to the 2D Anderson model at weak disorder p8| , 
where a qualitative change towards Poisson-like behav- 
ior for larger system sizes is observed, indicating a finite 
localization length of the eigenstates. The present size- 
independence of the LSD is compatible with multifractal 
and extended states. 

In conclusion, we have shown that the energy level 
statistics of TB Hamiltonians defined on the octagonal 
tiling with different boundary shapes is very well de- 
scribed by RMT. We can even see that our data fit the 
exact integrated LSD /goe(s) better than the integrated 
Wigner surmise iw( s )- This supports the applicability 
of RMT for such completely deterministic Hamiltonians. 
Although there is no randomness in these quasiperiodic 
models, one may view the absence of translational sym- 
metry as a sort of "topological disorder". We clarify 
previous statements p3)-^5| and show that the univer- 
sal LSD for irreducible blocks of a symmetric patch, or 
for patches without any symmetry is always Iqoe(s)- For 
patches with the full Z?s-symmetry, we find an integrated 
LSD which is a superposition of seven /goe(s)- Besides 
the octagonal tiling, we have also considered planar 10- 
and 12-fold quasiperiodic tilings and found analogous re- 
sults. In all these cases, we never find a critical I c (s), 
distinct from Igoe(s) and Ip{s), as observed at the An- 
derson MIT ^ . This is somewhat surprising since eigen- 
states in these quasiperiodic tilings are multifractal sim- 
ilarly to states at the MIT, and we could have expected 
that this is reflected in the LSD. Instead, we find that 
the LSD is similar to the LSD on the metallic side of the 
MIT. 
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FIG. 1. (a) Octagonal cluster of the Ammann-Bccnker 
tiling with exact eightfold rotational symmetry around the 
central vertex (x,y) — (0,0) as indicated by the solid and 
dashed lines. Shadings indicate successive inflation steps 
of the central octagon. The patch contains 833 vertices, 
(b) Square-shaped cut defined by < x < L, — -§■ < y < 
with L = 20 with 496 vertices, (c) Sinai billiard-shaped patch 
defined by < y < x < L and x 2 + y 2 > ^ with L = 22 with 
246 vertices. The gray edges correspond to the interior of the 
circular arc; edges crossing the arc have been deleted. 
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FIG. 2. Integrated LSD I(s) for (a) the largest 

D8-symmetric octagonal patch, crosses (x) indicate Ip(s), 

(7) 

diamonds (o) indicate 1^ (s); (b) the seven independent sub- 
spectra of the largest Ds-symmetric octagonal patch, circles 
(o) indicate 7w(s); (c) squared-shaped patches of different 
sizes, circles as in (b); (d) Sinai billiard-shaped patches of dif- 
ferent sizes, circles as in (b). Curves (b)-(d) have been shifted 
by multiples of 0.2 for clarity. 




FIG. 3. Small-s (left) and large-s (right) behavior of I(s) 
for one irreducible sector of Ds-symmetric octagonal patches 
of different sizes. The bold line corresponds to the largest 
patch. The three smallest and largest level spacings for each 
patch are denoted by triangles of different orientations. The 
circles (o) indicate Iw(s). Inset: blow-up of the data region 
enclosed by the rectangle, showing only data for the largest 
patch. Squares (□) indicate Jgoe(s). 
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